Correlated random fields in dielectric and spin glasses 
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PACS 64 . 70 . ph - Nonmetallic glasses 

PACS 75. 10. Nr - Spin-glass and other random models 

PACS 75 . 10 . Jm - Quantized spin models 

Abstract. - Both orientational glasses and dipolar glasses possess an intrinsic random field, 
coming from the volume difference between impurity and host ions. We show this suppresses 
the glass transition, causing instead a crossover to the low T phase. Moreover the random field 
is correlated with the inter-impurity interactions, and has a broad distribution. This leads to a 
peculiar variant of the Imry-Ma mechanism, with 'domains' of impurities oriented by a few frozen 
pairs. These domains are small: predictions of domain size are given for specific systems, and 
their possible experimental verification is outlined. In magnetic glasses in zero field the glass 
transition survives, because the random fields are disallowed by time-reversal symmetry; applying 
a magnetic field then generates random fields, and suppresses the spin glass transition. 



Most real solids are either amorphous or disordered 
crystals, and many of these show some sort of 'glassy' 
behavior at low temperatures [1] . These glasses have been 
categorized into different types, including spin [2J, struc- 
tural [3] , electron [3] , orientational [5J , and dipolar glasses 
[BJ. All of these show aging, rejuvenation and memory 
effects. However, thermodynamically there is a marked 
difference between systems like spin glasses (SGs), which 
show a clear phase transition between the paramagnetic 
and SG phases (with diverging length scale and conse- 
quent diverging nonlinear susceptibility [7]), and systems 
like orientational glasses (OGs) or electric dipole glasses 
(DG) where the transition to the low-T phase is smeared 

mm- 

At first glance, this difference seems paradoxical - in 
OGs, the rotational degree of freedom behaves like an 
Ising pscudospin 1/2, similar to that in DGs; and the 
electric- and strain-mediated interactions between these 
pseudospins have a similar form to those between the real 
spins in dipolar SGs. The quenched randomness then im- 
plies that OGs and DGs should show the same critical 
behavior as dipolar SGs, that of the short range Edwards- 
Anderson model [ID] . 

In the present paper we argue that the full answer to this 
question is found in an analysis of the effective Hamiltoni- 
ans for these systems - including the form of the random 
fields (RFs) that are generated - and in the constraints im- 



posed by symmetry. These RFs are broadly distributed, 
and correlated with the interactions; we find this leads to 
a peculiar disordering mechanism, intermediate between 
simple alignment along the RFs and Imry-Ma [TT] behav- 



ior. 



(i) Orientational &; Dipolar glasses: For typical 
OGs like KBr:CN, various approaches have been used to 
describe the interactions [T2TH5] . If one adds to the stan- 
dard strain-mediated interaction 14, 16 another term re- 
flecting the volume difference between host and impurity 
ions, an effective random field (RF) is found to be gener- 
ated by the orientational impurities themselves. The OG 
degrees of freedom can be treated formally in the same way 
as the pseudospin degrees of freedom in DG systems like 
Li:KCl, allowing one to analyse the interactions of both 
systems on the same footing [15] . 

One then finds |15j that over a wide range of temper- 
atures extending up to the putative glass transition tem- 
perature Tq, both OGs and DGs are described by the 
pseudospin effective Hamiltonian 



H. 



eff 
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We emphasize that both the Ising interaction U t 
the RFs {bj} are generated by the same defect-phonon 
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V dp = 
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expanded to second order in perturbation. Here Xj a 
represents the displacement operator at point j in di- 
rection a. One finds US] U zz = g zz T zz jR?^, and bj = 

Ei u °f = 9oz j2t '/,'r ''v with R y = r * ~ r,-, and where 

9zz ~ ll/pc 2 , and 3oz ~ (Vojo/pcf). Here p is the den- 
sity, c (c/) is an average (longitudinal) sound velocity, and 
T zz ,Tf z ~ 0(1) are complicated angular averages over 
the tensor components of the interaction. 

The typical parameters of the Hamiltonian (JTJ) are there- 
fore given by 



U 



P C o 



Volo 
2"' 



(3) 



and we note that the ratio U/b is independent of n. 

In principle, one should add to the hamiltonian ([I]) a 
tunneling term ^■j&j y which is crucial for the dynam- 
ics of the system. However, we are interested in the case 
where interactions dominate, ie., U > A; the tunneling 
term is then irrelevant when considering the phase dia- 
gram. For concentration n = 10 -3 , typical values are 
U ~ 1.5 K and b somewhat smaller. Typically we assume 
n <C 1 in what follows, but we will also compare with 
larger n. 

(ii) Correlated random fields: Consider now the 
RFs {bj}. These were previously argued to lead to a glassy 
state |13j . or to a state intermediate between a true spin 
glass and a pure random-field system |17) . Here we ar- 
gue that actually the RFs destroy long-range glassy or- 
der in the low-T phase [TSJ, which is then reached by 
a crossover and not a phase transition. Furthermore, 
since the RFs originate from the impurity lattice inter- 
actions themselves, they are correlated with the inter- 
actions (because U° z is correlated with U zz ) and have 
a broad distribution. When defects/impurities occupy 
nearest-neighbour sites, the field of one on the other is 
bo ~ Vo7o/(pcl) > b. 

Because Uij oc 1/i? 3 and is random, (Ufj) oc 1/R e and 
is effectively short range [TUJ[T5] . The situation is thus 
similar to that analysed by Imry-Ma [11], and its gen- 
eralisation to glasses by Fisher and Huse [15] . However 
there is an important difference - we find that in systems 
like dilute OGs, where the RFs are broadly distributed 
and are correlated with the interactions, one gets qualita- 
tively different results from systems where the RFs are 
normally distributed. The largest RF, occuring in the 
event when defects/impurities occupy nearest-neighbour 
sites, is &o- Since n -C 1, such events are rare; how- 
ever the average RF 6 av (the standard deviation of the 
distribution) is dominated by these rare events, and is 
given by 6 av ~ rj ^ ^/n/ (pc 2 ) » b (see Ref. [19], noting 



the trivial relations between high and low impurity densi- 
ties) . Now, since bo -C Uo (where Uq is the interaction be- 
tween nearest-neighbour defects/impurities), two regimes 
are possible, viz.: (a) bo <C Uqu and (b) Uon <C b <C Uo- 
Regime (a) is simple. All RFs are smaller than the typi- 
cal interaction. One can then follow the standard applica- 
tion of the Imry-Ma argument [11] to SGs [18] . The energy 
cost to flip a domain is related to the typical interaction, 
and therefore is ~ UL e , where L is the domain size, and 
the stiffness exponent 9 ~ 0.2 in 3D [T51I2D] . However, the 
effective RF is given by £» av = bo\/n = y/b~ob. Thus, the 
number of spins in a domain Nd, and the related dimen- 
sionless correlation length £ = ^n 1 / 3 
dependent, the latter being given by 



nV 3 are both n 
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In comparison to Gaussian-distributed RFs, the domain 
size here is smaller by a factor n ' 39 . 

Regime (b) (which usually applies in dilute dipolar 
glasses) is more complex, because the spread of RFs is 
larger than the typical interactions, so Imry-Ma argu- 
ments cannot be used directly. Consider a pair of nearest- 
neighbour impurities/defects 1,1' (3-impurity correlations 
can be neglected when n <C 1). Since Uo 3> U, the two 
impurities will mutually order, ie., lock together. If this 
locking is antiferromagnetic, then the net RF on the pair 
comes only from distant impurities - the internally gener- 



ated RF Uftf - U$f 



0, since U° z 



Uff for any pair i, j. 



However, for ferromagnetic ordering the internally gener- 
ated RF is ~ bo, and since Uon -C bo, this pair will be 
oriented along this internal RF, and hardly influenced by 
the interactions with other impurities, which have strength 
~ U = Uon\ In fact, this is the case for all ferromagnet- 
ically aligned impurity pairs satisfying bo/Rfj > U. In a 
sample of N impurities there are ~ Nbo/Uo such 'frozen 
pairs'. 

Consider now a sample containing N impurities/defects. 
The typical RF resulting from the frozen pairs is then 
~ \JboN /UoUon, while that from all other impurities is 
only ~ \/NbQn. Thus despite their scarcity, the RF from 
the frozen pairs dominates the total RF. Using now the 
standard argument which compares the energy gain from 
the RF and the interaction energy cost at domain bound- 
aries, we find for a dilute DG in regime (b) a dimensionless 
correlation length given by 



u 

lUb 



0.39 



(5) 



Thus, in this regime the correlation length is determined 
neither by the typical RF b n, nor by the average RF 
boy/n, but by the algebraic mean of U and b. This is 
because the relaxation of the RFs occurs in two stages. 
The rare pairs with > U relax locally to their frozen 
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state. The majority of the spins then order via the Imry- 
Ma mechanism, albeit with RFs dominated by their inter- 
action with the rare frozen pairs. Because U /b is scale- 
invariant Nd ~ (U/b) 1 - 17 , independent of n. The low 
power in the exponent is dictated by the unique mech- 
anism above, and results in "quasi-domains" , which in 
practise contain very few defects/impurities (see the ex- 
amples below). Moreover, in each such quasi-domain there 
are only ~ (U /b) a - 17 m 1 frozen pairs. In this limit one 
cannot of course assume that the random fields behave as 
in the large N limit; however, the mechanism described 
above gives an upper limit to the size of the domains, 
since it minimizes the effect of the RFs. A lower limit 
can be given by the assumption that the total RF from 
the frozen pairs is maximized, leading to Nd ~ (U/b). In 
practise these two limits are indistinguishable, suggesting 
that domains are formed in such a way that interactions 
with frozen pairs are maximized. 

Thus, in regime (b) we arrive at a picture of very 
small quasi-domains or 'clusters' of Nd pseudospins, of 
size £, each containing ~ O(l) frozen pseudospin, whose 
field then orients the rest of them. For n « 1 (regime 

(a) ) we find N d m (U/b) 2 - 34 , giving large clusters, but 
on strong dilution this crosses over to the much smaller 
N d ~ (U/b) 1 - 17 . We note that for all dilutions, the finite 
£ gives a finite nonlinear susceptibility X3 oa ^ 2_r ', mark- 
ing a crossover rather than a transition to the glassy state 
[3TJ[52], in agreement with experiments [5]. However, the 
decrease in Nd with n means that the peak in the nonlin- 
ear susceptibility should be much smaller for dilute OGs 
and DGs, in comparison with similar systems at high con- 
centrations. This is a central prediction of this Letter. 

In real systems, the ratio U /b varies widely between 
different OGs and DGs, and it is useful to look at some 
experimental examples. Consider first OH - impurities in 
KC1, where the electric and elastic impurity-impurity in- 
teractions are comparable, and U/b ~ 3 (see Ref. |23|). 
In Ref. [9] it was argued, based on the behavior of the 
non-linear dielectric permittivity of dilute KCl:OH, that 
there is no transition to a glass phase, but rather to a 
state analogous to a superparamagnet, with roughly 10 
impurities per domain. The scaling approach we use here, 
embodied in Eq. ([5]), is strictly applicable only for £ 3> 1, 
and even then only gives an order of magnitude. Thus 
our result here is consistent with the above experimen- 
tal picture of small 'superparamagnetic' domains for this 
system. Consider now systems where the spherical and 
dipolar species are of similar volume, so b/U -C 1. For 
example, b/U 1/20 in KBr : CN (see Ref. |Mj), which 
is an OG for n < 0.5. For n w 0.5 we obtain £ w 10 and 
Nd ~ (U/b) 2 - 34 w 1000. However, in the dilute regime 

(b) , because of the small exponent in Eq. ©, £ ~ 3, and 
Nd ~ 30 only, a large reduction! 

As mentioned above, the electric dipole interaction is 
typically quite small. For our purposes here, it renormal- 




E/n(*10 V/m) 



Fig. 1: The dimensionless correlation length f = £• n 1 ^ 3 of the 
glass order is plotted as function of the normalized external 
field E/n for KBr:CN (b/U = 1/20, d = 0.3 D) with n = 0.5 
(solid black line) and n — 10 -3 (dotted red line), against its 
value with zero effective RF (dashed blue line). The inset shows 
a log-log plot, emphasizing the n dependence of the deviation 
from power low behavior. 



izes the Ising interaction, so that 




(6) 



and allows the coupling to an electric field E. The latter 
allows a measurement of the effective RF as function of n, 
providing an experimental check to our theory. For E > 



6 



u 



0/2)-e 



^/b 2 cB + (E ■ df 



(7) 



where d is the dipole moment, b c g = boy/n in regime (a) 
and b c ff = VbU in regime (b). Thus, increasing E causes 
a crossover from impurity-dominated to field-dominated 
£ at E ■ d « b cS . In Fig. Q] we p_lot_ £ for KBr : CN, for 
n = 0.5 and for n = 10~ 3 < b/U. For E > b cS /d, 
one has £e oc E^Vts/s-e)^ However, the magnitude of 
£e at E = 0, as well as the region where its functional 
form deviates considerably from the above power low, are 
strongly n dependent. Their measurement as functions of 
n would measure the effective RF in the system as function 
of dilution, and thus provide a check to our theory. 

Extrinsic impurities also generate RFs with a broad dis- 
tribution, and can be analyzed along the same lines as the 
intrinsic impurities. Thus, at E = 0, as the system is pu- 
rified, the magnitude of £ and of the cusp in the nonlinear 
permittivity increase, but only to a value in accordance 
with Eqs. Q, ©. 

(iii) Spin glasses: We have seen that the random fields 
generated in OGs and DGs prevent a genuine phase tran- 
sition to a low-T glass state. Why does one then see a 
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glass transition in SGs? After all, one could certainly ex- 
pect similar random fields to be generated therein, via, 
eg., magnetoacoustic interactions. However, there is an 
essential difference between electric and magnetic Ising 
systems. If one neglects the volume term in Eq. @, then 
OGs share with SGs a symmetry under o z — > — a z . How- 
ever this is a symmetry under parity P in OGs (where 
the pseudospin variables a z are not real spins) , but under 
time-reversal T in SGs (where the a z are real spins, and 
couple to real magnetic fields). As a result, terms linear 
in a which emerge naturally in OGs, are not allowed in 
the absence of a magnetic field in SGs [35] ■ For example, 
in zero field the magnetoacoustic interaction is, to lowest 
order, bilinear in a: 



v sP = -EE 



a/3 



r)5 



a/3 
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dX 
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dx 
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where A^ lS is the spin-phonon interaction tensor. Thus, 
SGs are well described by the Edwards Anderson model 
with no RF, and have a well defined phase transition and 
diverging nonlinear susceptibility between the SG and PM 
phases [7]. 

The fact that it is the time-reversal nature of the Z2 
symmetry that prevents the emergence of random lon- 
gitudinal fields is best exemplified in anisotropic dipolar 
magnets. Recently it was shown that with the applica- 
tion of a transverse field H±, that breaks time reversal 
but by itself keeps the a z —¥ —a z symmetry, an effective 
RF h emerges via the intrinsic off-diagonal terms of the 
dipolar interaction [26]. Similar RFs would emerge from 
the magneto-acoustic interactions in the presence of H±. 
Thus, anisotropic dipolar magnetic systems in a trans- 
verse field are equivalent to dilute electric systems, albeit 
with a tunable effective RF. In the SG phase, these ran- 
dom fields result in a crossover rather than a quantum 
phase transition between the SG and PM phases as func- 
tion of transverse field [35]. Indeed, experiments on the 
LiHo x Yi_ x F4 have found [57] that the cusp in the nonlin- 
ear susceptibility is smeared, with cusps becoming smaller 
as temperature is decreased, and the applied magnetic 
field at the crossover, and thus the effective random field, 
are increased. The above mechanism for RFs in magnetic 
system also applies to the FM regime, making anisotropic 
dipolar magnets the first realization of the RFIM in a fer- 
romagnetic system [T91f2"FIf2~ni ■ Thus, measuring the RF 
in anisotropic dipolar magnets would be of much inter- 
est: this could be done, eg., in the same way as suggested 
above for OGs, by mapping b -O- h and E ■<-> Hn, where 
if 1 1 is the longitudinal magnetic field. 

All our considerations for the emergence of effective 
RFs, both in OGs and in magnetic systems with trans- 
verse magnetic field, are independent of the thermody- 
namic phase of the system, i.e. they also apply to the 
ferromagnetic/ferroelectric phases [15] . Thus, we argue 
that in disordered system with no time-reversal symmetry 



the Ising model is unstable to small perturbations, and 
is therefore not realizable (i.e., RFs will always emerge). 
In particular, one cannot realize the transverse field Ising 
model in easy-axis disordered magnetic systems by apply- 
ing a transverse magnetic field. 

In light of these conclusions for SGs, it is interesting 
to recall that models of 2-level systems (TLSs) interact- 
ing with phonons are commonly used to describe a wide 
variety of glasses at low T, not just OGs and DGs. Inso- 
far as the models of TLSs are applicable, and the volume 
term in the interaction with phonons causes a breaking of 
the cr z o — a z symmetry, then our theory predicts that in 
these systems RFs will emerge as well, and there will also 
be a crossover rather than a phase transition between the 
glass-ordered and disordered phases. To unambiguously 
test this, one could measure the non-linear dielectric sus- 
ceptibility in these systems - it should not diverge as a 
function of temperature. These considerations apply both 
for the glass transition, and for the low energy regime 
of interacting TLSs [30]. For a recent calculation of the 
strength of the TLS-TLS interactions in the two regimes in 
orientational glasses, and a discussion of its applicability 
in structural glasses, see ref. [3"3"] . 
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